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As recently shown, Fermi-LAT measurements of the diffuse gamma-ray emission from the Galaxy 
favor the presence of a smooth softening in the primary cosmic-ray spectrum with increasing 
Galactocentric distance. This result can be interpreted in terms of a spatial-dependent rigidity 
scaling of the diffusion coefficient. The DRAGON code was used to build a model based on 
such feature. That scenario correctly reproduces the latest Fermi-LAT results as well as local 
cosmic-ray measurements from PAMELA, AMS-02 and CREAM. Here we show that the model, 
if extrapolated at larger energies, grasps both the gamma-ray flux measured by MILAGRO at 
15 TeV and the H.E.S.S. data from the Galactic ridge, assuming that the cosmic-ray spectral 
hardening found by those experiments at about 250 GeV/n is present in the whole inner Galactic 
plane region. Moreover, we show as that model also predicts a neutrino emission which may 
account for a significant fraction, as well as for the correct spectral shape, of the astrophysical 
flux measured by IceCube above 25 TeV. 
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1. Introduction 

Cosmic-ray (CR) transport in the Galaxy and secondary y-ray diffuse emission arc generally 
modeled assuming that charged particle species undergo a statistically homogeneous and isotropic 
diffusion due to their scattering onto the Galactic magnetic field random fluctuations. This amounts 
to treat the diffusion coefficient as a spatially independent scalar function of the particle rigidity. 
Although this simplified approach does not stand on solid theoretical grounds - the random mag¬ 
netic field power is not expected to be homogenous and the presence of a regular field component 
should break the isotropy - nevertheless it provides a quite successful description of a wide number 
of experimental data. Indeed, CR spectra of primary and secondary CR species can accurately be 
reproduced in that conventional framework at least up to few hundred GeV/nucleon (see e.g. [1] 
for some recent updated results). 

It should be noted, however, that charged CR species provide a rather local probe of the Galac¬ 
tic environments. Locally detected secondary species, as the boron, which are generally used to 
tune the relevant propagation parameters, are produced within few kpc from the Earth, therefore 
they do not probe the conditions in the central region of the Galaxy. Instead, secondary y-rays, pro¬ 
duced by the interaction of CR nuclei (electrons) with the interstellar gas (radiation), offer a deeper 
insight of the Galaxy. Interestingly, although the conventional picture implemented with numerical 
codes like GALPROP [2, 3] reproduces the main features of the Galactic diffuse emission, several 
anomalies suggest that some main ingredients may be missing in that approach. 

We start noticing that conventional models cannot explain the large y-ray flux measured by the 
Milagro observatory from the inner GP region (\b\ < 2°, 30° < / < 65°) at 15 TeV median energy 
[4]. An optimized model which was proposed [5] to account for the EGRET GeV excess and 
came out to reproduce Milagro result as well, was subsequently excluded by Fermi-LAT [6]. This 
problem holds also for the most updated conventional models tuned to reproduce updated Fermi- 
LAT data [7] (see below). Another intriguing discrepancy was found between the predictions of 
conventional models and high-energy y-ray H.E.S.S. data for the Galactic ridge region (|Z| < 0.8°, 
\b\ < 0.3°) [8]. The diffuse spectrum measured by H.E.S.S. for 0.3 < E < 10 TeV in that region 
is well described by a power law with index 2.29 ± 0.07 ± 0.20. This is significantly harder than 
expected from 7To-decay if a CR spectral shape identical to that found in the solar neighborhood is 
assumed in the GC region. While in [8] the authors suggest such discrepancy to be originated by the 
proximity between the dense molecular gas with a CR source in that region, no compelling evidence 
supporting that interpretation has been provided and other explanations worth to be considered. 

It is important to point out that several troubles for conventional models arise also at lower 
energies: GALPROP based conventional models tuned to reproduce the diffuse y-ray over most 
of the sky (see e.g. [2, 3] and references therein) systematically underestimate the measured flux 
in the inner Galactic Plane (GP) region above few GeV. This problem, which was pointed-out by 
the Fermi collaboration itself [7], has been independently confirmed using updated Fermi data and 
point source catalogue (3FGL) [14]. 

Other independent hints of the inadequacy of the standard CR propagation scenario are com¬ 
ing also from the newly born high-energy neutrino astronomy. In the last few years the Ice- 
Cube collaboration claimed the detection of 37 neutrino events above 30 TeV corresponding to 
an excess with respect to the atmospheric background of 5.7 a [9]. The inferred flavor composi- 
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tion is compatible with a mixture of electronic, muonic and tauonic neutrino in equal amounts 
as expected if their origin were astrophysical. A recent analysis [10] of 641 selected events 
with vertices contained in the detector allowed the IceCube collaboration to lower the energy 
threshold and to measure the extraterrestrial diffuse neutrino spectrum to be 4> v , = 2.06 1 g 4 x 
10 -18 (£ v /10 5 GeV) 2 - 46±012 GeV -1 cm -2 sr -1 s -1 for 25 TeV < E v < 1.4 PeV. This spectrum 
is significantly softer than expected from most extragalactic sources (e.g. AGNs or y-ray bursts). 
This fact, together with a predominancy of down-going events and their concentration - though still 
not statistically significant - toward the GC region, suggests the presence of a significant Galactic 
component of the detected astrophysical neutrino flux [11]. Although the Galaxy is known to be a 
guaranteed source of a neutrinos originated from the hadronic interactions between the Galactic CR 
sea and the diffuse interstellar gas, conventional models predict an emission which is considerably 
smaller and softer than the flux measured by IceCube (see e.g. [13] and ref.s therein). 

In this contribution we discuss a possible common solution of all those anomalies. This is 
based on a scenario, recently presented in [14], featuring a radial dependence for both the rigidity 
scaling index 8 of the diffusion coefficient and the convective wind. Since the model assumes 8 
increasing with R and it is tuned to reproduce CR locally observed spectra, it predicts a hardening 
of CR propagated spectrum, hence of the secondary y-ray emission, in the inner Galaxy. The 
model also accounts for the CR proton and Helium spectral harding inferred from PAMELA [15], 
CREAM [16] and AMS-02 [17] at few hundred GeV/n. We will see as such a feature is actually 
required to get a consistent description of Fermi-LAT and upper-TeV y-ray data by Milagro and 
H.E.S.S. 

2. The KRA r setup 

We treat CR transport in the Galaxy by means of the DRAGON code [19, 20]. Similarly to other 
numerical CR transport codes (e.g. GALPROP [2, 3]), DRAGON solves numerically the diffusion 
equation for a given source term and null boundary conditions away from the Galactic disk. While 
DRAGON shares with GALPROP the public CR spallation cross sections and the interstellar gas and 
radiation distributions (other choices of those quantities are however possible in DRAGON), differ¬ 
ently from GALPROP public versions DRAGON allows, among several other innovative features, to 
adopt a spatial dependent diffusion coefficient and convective velocity. This feature has been used 
to implement the KRA 7 models discussed in this contribution (see also [14, 21]). Those models 
assume that the exponent 8 determining the rigidity dependence of the CR diffusion coefficient 

D(R,z,p)=D 0 (p/p Q ) m exp (z/z t ) (2.1) 

has the following Galactocentric radial dependence: 8(R) = AR + B where A = 0.035 kpc 1 
and B = 0.21 so that 8(R. : ) = 0.5. Concerning the vertical dependence of the diffusion co¬ 
efficient we assume Zt = 4 kpc (we checked that our results do not change significantly con¬ 
sidering larger values of Zt)- As we mentioned, the models also adopt a convective wind for 
R < 6.5 kpc with velocity Vcizjz (z is the distance from the GP) vanishing at z = 0 and growing as 
dVjdz = 100 km s -1 kpc -1 . Those model parameters have chosen to consistently reproduce the 
slope and angular distribution of the diffuse y-ray emission of the Galaxy measured by Fermi-LAT 
as well as the main local CR observables (proton and Helium spectra and the B/C ratio) [14]. 
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Protons and Helium 



Figure 1: The computed proton and Helium (in the insert) spectra at the solar circle are compared with 
a representative set of experimental data. All spectra account for solar modulation with a potential <t> = 


0.5 GV. 


Since Fermi-LAT measured the spectrum of the diffuse Galactic emission up to few hundred 
GeV, some further hypothesis has to be done about the spectrum of primary CR above the TeV. 
In particular, concerning the p and He spectral hardening inferred from PAMELA [15] - recently 
confirmed by AMS-02 [17] - and CREAM [16] data above ~ 250 GeV/n, we consider two cases: 
1) KRAy no hard, model: the hardening is assumed to be a local effect - e.g. originated by 
nearby supernova remnants - which averages out on large scales hence introducing no feature in 
the Galactic CR population. 2) KRAy model: in this case the hardening is taken to be representative 
of the whole Galactic population. In both cases we assume that above 250 GeV/n the CR source 
spectra extend steadily up to an exponential cutoff varying in the range 5 < E cut /Z < 50 PeV to 
roughly agree with KASCADE [22] and KASCADE-Grande data [23] (see fig. 1). 


3. The diffuse 7 -ray spectrum above the TeV 

It was shown in [14], that both the KRAy no hard, and KRAy models provide good fits of the 
7 -ray diffuse emission measured by Fermi-LAT all over the sky, in particular towards the inner GP 
region. In [21] we extended that computation above the TeV. Starting from the propagated proton 
and Helium distributions we computed the hadronic emission integrating the expression of the 7 - 
ray emissivity [26] along the line-of-sight. For the gas components we used the distributions taken 
from the latest public version of the GALPROP package [2, 7]. 

In Fig. 2 we compare our predictions for the KRAy and KRAy no hard, with Fermi-LAT 
and Milagro diffuse emission data from the same inner GP region. In the same plot we also re¬ 
port DRAGON predictions for a representative conventional model (KRA) - again with and without 
hardening - sharing the same main properties of the KRAy model at the solar circle (;. e. 8 = 0.5 
and dV c /dz = 0 in the whole Galaxy). 
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30° < III < 65° Ibl < 2° 



x 10 12 Milagro longitude profile 



Figure 2: The diffuse emission y-ray spectrum (left panel) from the inner Galactic plane (\b\ < 2°, 30° < l < 
65°) computed for the reference models considered in this Letter is compared with Fermi-LAT and Milagro 
data. The expected sensitivity of HAWC and CTA are reported. The spectral components are shown for the 
KRAy model only. In the right panel the longitude profile at 15 TeV for the same models is compared with 
HAWC data. 


The reader can see as conventional models, which under-predict Fermi-LAT data in agreement 
with [14], more seriously fail reproducing Milagro result. While the presence of hardening in the 
whole CR Galactic population clearly enhances the y-ray flux above few hundred GeV, this is not 
enough to account for the flux measured by Milagro from the inner GP at 15 TeV. The KRAy setup, 
instead, is more successful - especially if the CR hardening is assumed. In the right panel of Fig.2 
we also compare our results with the longitude profile at 15 TeV measured by Milagro. Although 
the data do not allow any strong claim, we notice a better fit of the KRAy model, although this 
should only be taken as a hint. 

Exploring a region not covered by Milagro, the Fl.E.S.S. observation of the Galactic ridge 
offers a further test of our model. Interestingly, also for that region the model provides a good 
fit of Fermi-LAT and Fl.E.S.S. results (see Fig. 3). We also cross-check adopting a more realistic 
gas distribution [24] in the inner Galaxy, and we rescale the models by a factor of 0.3 to minimize 
the x 2 against the data: This factor is justified by the smaller value of the conversion factor (Vco) 
between the FE column density and the CO line brightness temperature in the central region [24]. 

4. Implications for the Galactic neutrino emission 

The presence of a hard CR spectrum in the inner Galactic plane region as predicted by the 
KRAy setup has relevant consequences also for the diffuse neutrino emission produced by CR 
hadron scattering onto the interstellar gas. In order to estimate such an effect, we computed the 
v e and v u production spectra similarly to what done for y-rays. For both flavors we used the 
emissivities provided in [26] for projectile energies below ~ 500 TeV, while we adopted the ones 
provided in [27] above that range. We also accounted for neutrino oscillations which redistribute 
the composition almost equally among all the three flavors. We only considered proton and Helium 
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III < 0.8° Ibl < 0.3° 



Figure 3: The computed 7 -ray diffuse emission from the Galactic ridge region is compared with Fermi-LAT 
and H.E.S.S. data. For each model the spectrum normalization has been varied to minimize the x~ against 
the data. The spectral components are shown for the KRAy model only. 


CRs/gas - as done for 7 -rays - since heavier nuclear species give a negligible contribution in the 
energy range we cover in this work. 

In Fig. 4 we show the total neutrino spectrum (all flavors, including antiparticles) computed 
for the KRAy (blue lines) in the presence of a CR spectral hardening in the whole Galaxy. In order 
to account for the uncertainties on the high energy tail of the primary CR spectra we considered 
two values of the CR source spectrum cutoff £ cut = 5 PeV (dashed lines) and 50 PeV (solid lines) 
roughly bracketing KASCADE and KASCADE-Grande data (see Fig. 1). For comparison we also 
show our predictions for the convetional KRA model for the same values of E CM (red lines). In all 
cases we assume the CR hardening at ~ 250 GV to be present in the whole Galaxy. The reader 
can see from that figure that, while the KRA setup predicts an undetectable flux - in agreement 
with previous results for other conventional models [19] - the KRAy model is instead able to 
explain a significant fraction of the flux measured by IceCube in terms of Galactic diffuse emission. 
Interestingly, for that model, the flux spectrum below 100 TeV - where the Galactic contribution is 
expected to be larger - almost follows a power law spectrum with an index which is very close to 
that measured by IceCube. 

Setting a threshold energy at 25 TeV, with a cosmic-ray source cutoff E cul = 50 PeV we 
obtain an expected event rate representing 40% of the complete sample of 37 events reported, well 
above the expected rate of atmospheric muons and atmospheric neutrinos after the veto conditions 
imposed by the IceCube collaboration. On the other hand, considering a threshold energy of 60 
TeV the expected event rate is about 30% of the corresponding measured sample of 20 events. 
Decreasing the E cut to 5 PeV the expected event rates become 20% (above 25 TeV) and 10% 
(above 60 TeV) of the measured IceCube samples. This result could soon be confirmed for the 
southern hemisphere emission by ANTARES experiment [28, 29] selecting a low energy sample 
of events (< 60 TeV) and smaller portion of the sky, and accurately verified when IceCube and the 
incoming KM3NeT [30] obtain a large catalog of astrophysical neutrino events. 
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Figure 4: The full-sky neutrino spectrum computed for the KRAy and KRA models (both with global CR 
hardening) and two choices of the CR high-energy cutoff is compared with IceCube result [10]. 


Concerning the angular distribution of IceCube events, we notice that - although IceCube col¬ 
laboration did not report any statistically significant anisotropy yet - a recent analysis [31] showed 
as a Galactic component producing up to 50% of the astrophysical neutrino flux measured by Ice- 
Cube is compatible with the angular distribution of high-energy events detected by that experiment. 
Therefore, our model is consistent with those results. Interestingly, the IceCube collaboration just 
published a new analysis [32] pointing to different flux normalizations and spectral indexes in the 
North and South hemisphere which may be interpreted with the presence of a component from the 
inner Galaxy as suggested in the above. 


5. Conclusions 

Conventional CR transport models - assuming uniform diffusion/advection properties - pro¬ 
vide unsatisfactory description of the diffuse Galactic y-ray diffuse emission measured by Fermi- 
LAT toward the inner Galactic plane. Moreover these models cannot reproduce measurements 
taken in the same region above the TeV. We showed that a recently proposed scenario assuming 
a proper radial dependence for both the rigidity scaling index 8 of the diffusion coefficient and 
the convective wind, which allows to reproduce Fermi-LAT results [14], also matches Milagro 
observations at 15 GeV. Using a detailed gas model for the Galactic center regions, we showed 
as that model also consistently reproduces Fermi-LAT and H.E.S.S. measurements in the Galac¬ 
tic ridge region. Those results requires the CR spectral hardening found by PAMELA, AMS-02 
and CREAM experiments at about 250 GV to be not a local effect but to be present also in inner 
Galactic plane region. We also considered the implications of our scenario for the neutrino diffuse 
emission of the Galaxy, Interestingly, we found that depending on the details of the CR spectra in 
the knee region, our model predicts that a fraction between 10 % and 40% of the signal measured 
by IceCube may be originated by CR scattering with the Galactic interstellar medium. 
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